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HIGHLIGHTS 


•  Extended  accelerated  aging  tests  on  lithium-ion  batteries  including  storage  and  cycling. 

•  Detailed  analysis  of  temperature  and  voltage  dependencies  of  calendar  cell  aging. 

•  Detailed  analysis  of  influence  of  cycle  depth  and  SOC  range  on  cycle  aging. 

•  Lifetime  prediction  based  on  identified  main  aging  phenomena. 

•  Comprehensive  verification  measurements  for  aging  model. 
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Knowledge  on  lithium-ion  battery  aging  and  lifetime  estimation  is  a  fundamental  aspect  for  successful 
market  introduction  in  high-priced  goods  like  electric  mobility.  This  paper  illustrates  the  parameteri¬ 
zation  of  a  holistic  aging  model  from  accelerated  aging  tests.  More  than  60  cells  of  the  same  type  are 
tested  to  analyze  different  impact  factors.  In  calendar  aging  tests  three  temperatures  and  various  SOC  are 
applied  to  the  batteries.  For  cycle  aging  tests  especially  different  cycle  depths  and  mean  SOC  are  taken 
into  account.  Capacity  loss  and  resistance  increase  are  monitored  as  functions  of  time  and  charge 
throughput  during  the  tests.  From  these  data  physical  based  functions  are  obtained,  giving  a  mathe¬ 
matical  description  of  aging.  To  calculate  the  stress  factors  like  temperature  or  voltage,  an  impedance 
based  electric-thermal  model  is  coupled  to  the  aging  model.  The  model  accepts  power  and  current 
profiles  as  input,  furthermore  an  ambient  air  temperature  profile  can  be  applied.  Various  drive  cycles  and 
battery  management  strategies  can  be  tested  and  optimized  using  the  lifetime  prognosis  of  this  tool. 
With  the  validation  based  on  different  realistic  driving  profiles  and  temperatures,  a  robust  foundation  is 
provided. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium-ion  batteries  are  a  key  technology  for  current  and 
future  energy  storage,  whether  they  are  used  for  mobile  or  sta¬ 
tionary  application.  As  the  batteries’  portion  of  cost  is  quite  high  for 
many  applications,  battery  lifetime  is  a  critical  point  for  profit¬ 
ability.  However,  performing  real  life  aging  tests  for  every  single 
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application  is  an  expensive  and  time  consuming  process  which 
cannot  be  done  in  every  case.  Lifetime  prediction  using  aging 
models  can  overcome  this  challenge  as  they  have  to  be  done  only 
once  per  cell  type. 

Different  approaches  are  used  for  modeling  of  lithium-ion  bat¬ 
teries.  The  most  realistic  models  are  physico-chemical  ones  which 
simulate  mostly  single  aging  effects.  However,  these  models  are 
very  slow  and  complex  to  parameterize.  Therefore  they  are  not 
suitable  for  aging  prediction  on  a  long  time  scale.  In  contrast 
neuronal  network  models  have  a  very  high  computation  speed,  but 
they  are  usually  not  able  to  extrapolate  aging  for  conditions  that 
were  not  included  in  the  learning  set.  Using  physical  based  func¬ 
tions  fitted  to  accelerated  aging  tests  is  a  compromise  between  fast 
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models  on  the  one  and  realistic  models  on  the  other  hand.  These 
mathematical  functions  allow  to  reproduce  the  factors  influencing 
aging,  such  as  temperature  or  storage  voltage,  and  extrapolating 
from  a  fixed  set  of  tests  to  a  wide  range  of  applications. 

For  an  aging  prediction  a  lot  of  different  aging  factors  have  to  be 
accounted.  These  factors  are  e.g.  temperature,  storage  voltage  and 
time  for  calendar  aging  and  in  addition  cycle  depth,  SOC  range, 
current  rate  and  charge  throughput  for  cycle  aging  [1,2],  Some 
studies  have  already  published  simulations  dealing  with  a  few  of 
these  factors,  either  calendar  life  [3—6]  or  cycle  life  [7],  Completely 
parameterized  models  including  all  of  these  factors  and  consid¬ 
ering  both  calendar  and  cycle  life  are  still  missing. 

This  paper  shows  an  accelerated  aging  test  set  including  more 
than  sixty  cells  of  a  commercial  high-energy  18650  system  with 
NMC  cathode  material.  Both  capacity  loss  as  well  as  resistance  in¬ 
crease  are  addressed.  Calendar  and  cycle  aging  are  considered 
separately.  A  holistic  aging  model  including  an  impedance  based 
electric-thermal  model  is  presented.  Finally  the  comparison  be¬ 
tween  simulation  results  and  verification  measurements  at  nine 
different  conditions  will  be  shown. 

2.  Experimental 

The  tested  battery  was  the  Sanyo  UR18650E,  an  18650  round  cell 
which  is  manufacturer  rated  with  2.05  Ah  minimum  and  2.15  Ah 
typically.  The  cathode  active  material  is  Li(NiMnCo)C>2  (NMC)  and 
the  anode  consists  of  graphite.  It  is  a  high  energy  cell  with  a 
maximum  discharge  rate  of  3C  and  165  Wh  kg  1  energy  density. 
Voltage  limits  are  2.5  V  for  discharging  and  4.2  V  for  charging  with  a 
proposed  end  of  discharge  voltage  of  3.0  V  and  end  of  charge 
voltage  of  4.1  V. 

An  OCV  curve  of  the  cell  is  shown  in  Fig.  1.  A  detailed  analysis  of 
the  cell,  the  aging  tests  and  the  occurring  processes  can  be  found 
from  Ecker  et  al.  [8], 

2.1.  Calendar  aging  matrix 

A  list  of  all  tested  calendar  aging  test  conditions  can  be  found  in 
Fig.  2.  As  the  temperature  dependency  with  an  Arrhenius  equation 
is  already  described  by  some  authors  [1,9],  only  three  different 
temperatures  are  tested  to  verify  this  dependency.  More  attention 
was  given  to  the  voltage  dependency,  which  is  highly  resolved  with 
10  different  conditions.  Each  combination  of  temperature  and 
voltage  in  the  calendar  aging  test  consists  of  three  cells  to  get  some 
statistics.  All  cells  were  held  at  a  constant  voltage  (float  conditions) 
during  the  tests.  Every  seven  weeks  a  measurement  of  the  cells 
capacity  and  inner  resistance  was  performed. 

2.2.  Cycle  aging  matrix 
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Fig.  2.  Test  matrix  of  calendar  aging  tests. 


can  be  found  in  Fig.  3.  All  cycling  tests  were  done  at  a  current  rate  of 
1 C  and  a  cell  temperature  of  35  °C.  Temperature  was  logged  using  a 
sensor  mounted  on  the  surface  of  the  cell.  The  mean  SOC  was  set  by 
a  constant  voltage  charge  using  the  OCV  curve  of  the  cell.  Cycling 
was  done  Ah  based  around  this  point,  with  a  reset  of  the  mean  SOC 
every  100  cycles.  If  a  test  had  0  or  100%  SOC  as  limit,  these  points 
were  used  to  set  the  cycle  range.  Checkups  measuring  the  capacity 
and  the  inner  resistance  of  the  cell  were  made  every  3  weeks. 

2.3.  Checkups 

All  capacity  and  resistance  measurements  were  done  at  35  °C. 
The  capacity  measurement  starts  with  a  standard  charge  which  is  a 
1 C  charge  up  to  4.2  V  followed  by  a  CV  charge  until  current  was 
below  C/50.  After  that  the  capacity  was  measured  during  the  1C 
discharge  down  to  2.5  V.  The  inner  resistance  was  measured  at  SOC 
steps  of  10%,  starting  from  90%  SOC  down  to  10%  SOC.  Every  step  was 
reached  Ah  based  by  discharging  1/10  of  nominal  capacity  starting 
from  a  completely  charged  battery  in  the  first  step.  At  each  step  a 
pulse  power  characterization  profile  (PPCP)  is  applied  to  the  battery. 

The  PPCP  consists  of  an  18  s  2C  discharge  followed  by  a  40  s  rest 
period.  After  that  a  10  s  1C  charge  is  applied,  again  followed  by  a 
40  s  rest  period.  A  voltage  response  to  this  PPCP  is  shown  in  Fig.  4. 
From  this  profile  various  inner  resistances  are  calculated,  a  2, 10  and 
17  s  resistance  for  the  discharge  and  a  2  and  10  s  resistance  for  the 
charge.  For  the  aging  calculation  the  10  s  discharge  resistance  at 
50%  SOC  is  used.  This  resistance  is  calculated  as  the  difference  be¬ 
tween  the  voltage  before  the  discharge  pulse  and  10  s  after  the 
beginning  of  the  discharge  pulse  divided  by  the  current.  All  other 
resistances  are  calculated  in  a  similar  way. 

3.  Calendar  aging 

In  the  calendar  aging  tests,  cells  were  stored  at  different  tem¬ 
peratures  and  voltages.  Each  test  condition  was  performed  with  3 


The  main  focus  of  the  cycle  aging  tests  was  the  influence  of  cycle 
depth  and  mean  SOC.  A  diagram  with  all  cycle  aging  test  conditions 
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Fig.  1.  OCV  i 


1  for  a  new  cell  (Sanyo  UR18650E),  measured  at  35  °C. 


Fig.  3.  Test  matrix  of  cycle  life  tests  performed  on  a  2.05 
were  done  at  1C  and  a  cell  temperature  of  35  °C. 
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Fig.  4.  Voltage  reaction  of  a  new  cell  to  a  PPCP.  The  two  voltages  used  for  calculating 
the  10  s  discharge  resistance  are  marked  red.  (For  interpretation  of  the  references  to 
colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


cells  to  demonstrate  the  reproducibility  of  the  experiment.  The 
results  show  a  very  similar  aging  for  cells  tested  under  the  same 
conditions,  the  measured  capacities  show  an  especially  great  uni¬ 
formity  [8], 

Only  the  tests  at  100%  had  a  difference  of  more  than  4%  points 
between  the  best  and  worst  cell.  As  these  cells  also  had  a  very 
strong  aging,  they  are  excluded  from  further  analysis.  With  a  stor¬ 
age  voltage  of  4.162  V  they  were  above  the  recommended  end  of 
charge  voltage.  This  might  lead  to  additional  aging  effect  which 
cannot  be  scaled  down  to  lower  voltages. 

For  all  tests  a  mean  capacity  loss  and  resistance  increase  has 
been  calculated  for  every  checkup.  An  error  bar  plot  of  both  ca¬ 
pacity  and  resistance  in  the  50  °C  tests  can  be  found  in  Fig.  5.  The 
cells  suffer  increasing  capacity  loss  and  resistance  increase  with 
higher  storage  SOC  and  thereby  higher  voltage.  Also  the  tests  not 
displayed  here  with  different  temperatures  show  the  expected 
increased  aging  with  higher  temperatures.  A  detailed  discussion 
can  be  found  in  Ref.  [8], 

3.1.  Fit  method 

To  obtain  a  lifetime  model  out  of  the  accelerated  aging  tests,  a 
mathematical  description  of  the  aging  processes  is  needed.  The 
analysis  was  split  into  calendar  and  cycle  aging.  In  the  calendar 
aging  tests,  temperature  and  storage  voltage  were  varied;  a  func¬ 
tion  should  include  these  factors. 

From  Fig.  5  it  can  be  seen,  that  different  test  conditions  have 
unequal  checkup  numbers.  This  is  due  to  the  removal  of  cells  which 
went  below  75%  of  initial  capacity  or  a  later  start  after  the  first  tests 
were  done.  For  an  aging  function  being  directly  fitted  to  all 
checkups  this  would  have  lead  to  variable  impacts  of  different  test 
conditions  on  the  result.  The  influence  of  tests  with  moderate  aging 
would  have  been  overestimated  while  new  started  tests  or  tests 
with  higher  aging  and  therefore  fewer  checkups  would  have  been 
underestimated. 

To  avoid  this  a  two  step  fitting  was  done.  In  the  first  step,  every 
single  aging  condition  was  fitted  separately.  For  this,  different 
fitting  functions  were  tested.  The  best  choice  were  functions,  that 
had  only  one  aging  factor  as  fitting  parameter.  Afterward  the  de¬ 
pendency  of  this  aging  factor  on  temperature  and  voltage  was 
analyzed. 

3.2.  Time  dependency 

For  a  cell  with  a  carbon  based  anode,  like  the  one  used  here,  it  is 
widely  accepted  in  literature  [10,11]  that  the  dominant  calendar 
aging  effect  is  the  formation  of  a  so  called  solid  electrolyte 


time  /d 


Fig.  5.  a)  Normalized  capacity  over  time  and  b)  normalized  resistance  over  time  for 
calendar  aging  tests  at  50  °C.  For  each  SOC  the  average  on  three  cells  under  tests  is 
shown.  The  capacities  for  cells  under  same  conditions  have  a  great  uniformity,  while 
resistance  increase  varies  a  bit. 


interphase  (SEI).  The  SE1  is  built  up  of  decomposition  products  of 
the  electrolyte,  consuming  lithium  during  formation  and  increasing 
resistance  through  the  growing  layer  thickness  [12],  Although  there 
are  different  theories  on  SEI  formation  [  13,14],  both  lead  to  a  square 
root  of  time  dependency.  For  some  already  published  aging  tests 
[9,13]  this  seems  to  be  approved.  Nevertheless  the  capacity 
decrease  measured  at  this  cell  started  with  lower  aging  and  seemed 
to  have  a  share  of  linear  aging  which  was  also  observed  in  Ref.  [15], 
Therefore  a  pure  linear  and  a  pure  square  root  function  are  used 
together  with  a  superposition  of  both  and  a  function  with  t0,75  as  a 
medium  function.  The  fits  are  done  for  both  capacity  C  and  resis¬ 
tance  R,  using  n  and  a  (or  ai  and  a2)  as  fitting  variables. 


C/R  = 

n+a-t 

(1) 

C/R  = 

n+a-Vt 

(2) 

C/R  = 

n+a-lt+a2-Vt 

(3) 

C/R  = 

n  +  a-t075 

(4) 

For  all  fit  functions  the  parameter  n  should  be  equal  to  1  as  all 
data  is  normalized  to  initial  values.  Anyway  this  parameter  is  set 
free  for  the  fit  to  compensate  deviations  of  the  measured  data.  For 
further  analysis  only  the  aging  factors  a  are  used. 

An  example  of  the  fit  functions  can  be  found  in  Fig.  6.  The  mean 
capacity  of  cells  stored  at  50  °C  and  50%  SOC  is  plotted  together 
with  the  fit  functions,  showing  a  bad  congruence  of  single  linear  (1 ) 
and  square  root  functions  (2).  The  superposed  linear  and  square 
root  function  (3)  and  the  function  with  t075  (4)  are  similar  and 
match  the  data  well.  The  fit  parameters  for  50  °C  and  50%  can  be 
found  in  Table  1  together  with  an  average  R2  for  all  tests.  In  this 
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Fig.  6.  Mean  capacity  for  the  cells  stored  at  50  °C  and  50%  SOC  over  time  together  with 
4  different  fit  functions. 

table  and  all  following  equations  the  fitting  parameters  are  given 
without  units.  Instead  the  units  of  the  input  values  are  specified  to 
allow  recalculation.  For  Table  1  all  times  t  has  to  be  used  in  days. 

Best  results  on  the  single  tests  were  obtained  using  the  super¬ 
posed  function  (3),  but  the  aging  factors  a\  and  a2  do  not  have  a 
clear  trend  over  voltage  (see  Fig.  7).  This  might  be  due  to  scattering 
of  the  data  which  prohibits  an  exact  division  into  linear  and  square 
root  functions. 

Therefore  the  function  with  t0  75  is  used  for  further  analysis  of 
the  aging  parameters.  This  function  shows  a  similar  shape  in  the 
required  range  and  has  only  one  aging  factor  which  needs  to  be 
evaluated.  The  average  R2  value  is  not  as  good  as  the  one  from  the 
superposed  function,  but  better  than  the  single  linear  or  square  root 
functions. 


3.3.  Voltage  dependency 

A  plot  of  the  aging  factor  over  voltage  for  all  tests  performed  at 
50  °C  both  for  resistance  and  capacity  can  be  found  in  Fig.  8.  The 
aging  factor  show  a  linear  trend  considering  the  whole  voltage 
range.  Therefore  the  function 

av(V)  =  avV  +  a2  (5) 

is  chosen  here  to  describe  the  dependency  of  aging  factor  a  from 
voltage  V,  using  the  fitting  parameters  cq  and  a2. 

In  contrast  other  papers  found  an  exponential  dependency  [2] 
that  is  clearly  not  given  here.  Other  work  also  done  on  these  data 
sets  [8]  proposed  a  correlation  between  aging  and  graphite  stages 
in  the  anode  material,  dividing  the  voltage  dependency  in  single 
steps  with  nearly  equal  aging. 

In  this  case  for  the  aging  model  the  small  differences  between 
a  stepwise  function  and  a  simple  linear  dependency  can  be 


3.3  3.4  3.5  3.6  3.7  3.8  3.9  4  4.1 

voltage  /V 


Fig.  7.  Voltage  dependency  of  the  aging  factors  a)  a,  and  b)  a2  over  voltage  from  the 
superposed  linear  and  square  root  function  for  tests  performed  at  50  °C. 


neglected.  For  most  applications  voltage  and  SOC  are  changing 
during  the  simulated  profile,  sweeping  over  areas  which  underes¬ 
timate  or  overestimate  the  aging.  In  total  these  effects  should 
compensate  over  the  full  range.  Therefore  a  linear  regression  was 
used  to  describe  the  voltage  dependency  in  this  model. 

In  the  future  a  better  knowledge  on  the  aging  behavior  at 
different  graphite  stages  could  lead  to  a  simplified  parameteriza¬ 
tion  of  a  stepwise  defined  function  with  a  reduced  number  of 
needed  tests. 

3.4.  Temperature  dependency 

For  temperature  dependency  only  three  different  temperatures 
were  tested.  A  lot  of  literature  already  described  the  Arrhenius 
equation 


aT(T)  =  arexp(-^j  (6) 

as  a  good  matching  function  for  temperature  dependency.  Fa  is 
thereby  the  activation  energy  of  a  reaction  happening  at  a  tem¬ 
perature  T,  R  is  the  gas  constant.  To  test  this  equation  with  our  data, 
a  plot  of  ln(a)  over  inverse  temperature  is  shown  in  Fig.  9.  The  three 
temperatures  were  35  °C,  40  °C  and  50  °C.  For  both  capacity  and 
resistance  the  logarithm  of  aging  factors  give  a  linear  dependency, 
confirming  the  applicability  of  the  Arrhenius  equation  for  these 


Fit  parameter  for  capacity  and  resistance  from  the  50  °C  and  50%  SOC  test.  Four  different  fit  functions  with  three  to  four  free  parameters  were  used.  0R2  is  the  mean  R2  value  of 
all  fitted  tests.  All  fitting  values  are  given  without  units,  time  values  must  be  used  in  days. 


Capacity  Resistance  Capacity  Resistance  Capacity  Resistance  Capacity  Resistance 


0.982  1.045  1.020  0.981 

0.00035  0.00055  0.00862  0.01402 


R2  0.978  0.949  0.975  0.988 

0R2  0.974  0.952  0.960  0.965 


1.003  0.995 

0.00018  0.00015 

0.00427  0.01037 

0.997  0.994 

0.997  0.987 


0.999  1.017 

0.00173  0.00278 
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«T, scaled^)  =  Td).  £  (10) 

The  dependencies  of  the  two  single  fits  were  combined  and 
again  scaled. 

a(T,  V)  =  (r)JV--'ed(U)  (12) 

The  resulting  combined  equation  is 


voltage  /V 


Fig.  8.  Voltage  dependency  of  the  aging  factor  a  over  voltage  from  C  =  n  -  a  t0  75  for 
tests  performed  at  50  °C  together  with  a  linear  regression.  Although  the  resistance  has 
higher  aging  factors,  capacity  aging  is  restricting  lifetime  here  as  end  of  life  is  often 
defined  as  20%  of  capacity  loss  or  100%  resistance  increase. 


(T  _  av{V)-aT(T)-a0 
1 ’  ’  av(3.699V)-aT(50°C)' 


(13) 


The  combined  aging  factors  for  capacity  and  resistance  are 
acap  (T,V)  =  (7.543-V  —  23.75)-106-e_6?s  (14) 


cells.  Results  of  the  linear  interpolation  are  activation  energies  of 
58.0  kj  mol-1  I<-1  for  the  capacity  and  49.8  kj  mol-1  K  1  for  the 
resistance.  Therefore  a  change  from  20  to  30  °C  accelerates  calendar 
aging  by  a  factor  of  2.2  for  capacity  loss  and  2.0  for  resistance 
increase. 


3.5.  Fit  function  calendar  aging 


For  a  mathematical  model  of  calendar  aging  the  dependencies 
on  voltage  and  temperature  need  to  be  combined.  Here  an  aging 
factor  a  is  used  to  describe  the  aging  rate  during  a  test  of  time  t. 

C  =  t  —  acap  •  t0  75  (7) 

R  =  l  +  «res-ta75  (8) 

The  test  at  50  °C  and  50%  SOC  (3.699  V)  is  included  both  in 
voltage  (cfv.  see  Eq.  (5))  and  temperature  (or,  see  Eq.  (6))  fit  and 
therefore  chosen  as  an  intersection  for  both  functions.  As  the  two 
fitting  functions  do  not  exactly  match  at  this  point,  the  resulting 
mean  value  was  calculated. 


a(50°C,  3.699  V) 


av(3.699V)  +  aT(50°C)  _ 


(9) 


Both  functions  were  then  scaled  to  match  this  mean  value. 


inverse  temperature  /K  1  x  10 


Fig.  9.  Arrhenius  plot  of  aging  factor  a  for  both  capacity  and  resistance  over  inverse 
temperature.  The  shown  temperatures  are  35  °C,  40  °C  and  50  °C.  Also  the  linear 
regression  for  temperature  dependency  is  shown. 


ares(T,V)  =  (5.270-V-  16.32)-105-e“5TS  (15) 

for  tests  done  at  a  voltage  V  in  volts  and  an  absolute  temperature  T 
in  kelvins.  These  equations  allow  one  to  calculate  the  calendar 
aging  for  arbitrary  conditions  of  voltage  and  temperature.  They  are 
the  basis  for  the  cycle  aging  analysis  in  the  next  chapter. 

4.  Cycle  aging 

Cycling  a  battery  leads  to  additional  aging  due  to  processes 
which  do  not  occur  during  calendar  aging.  During  intercalation  and 
de-intercalation  the  material  experiences  a  volume  change  which  is 
a  stress  factor  for  the  battery  system.  Results  can  be  a  crack-and- 
repair  of  the  SEI  which  consumes  lithium  and  increase  the  inner 
resistance  or  a  contact  loss  of  active  material  particles. 

22  cycle  aging  tests  varying  cycle  depth  and  average  SOC  were 
performed.  All  cycle  aging  tests  were  done  at  a  temperature  of  35  °C 
and  a  current  of  1 C.  This  would  result  in  12  equivalent  full  cycles  per 
day,  but  real  values  are  lower  due  to  time  for  checkups  and  other 
unplanned  rest  periods. 

During  the  time  of  cycling  calendar  aging  also  occurs.  To  get  an 
analysis  of  the  ’pure’  cycle  aging,  all  test  data  had  to  be  adjusted  by 
the  calendar  aging.  The  calendar  capacity  loss  and  resistance  in¬ 
crease  were  calculated  using  the  functions  obtained  from  the  tests 
discussed  before  and  then  added/subtracted  from  the  measured 
capacities  and  resistances. 

Ccyc  =  Cmeas  +  «cap(Vmeas,  Tme as)  •  t0J5  (16) 

Rcyc  =  Rmeas  ~  «res(Vmeas,  Was)  •  t°J5  (17) 

To  calculate  the  calendar  aging  of  a  cycle  aging  test,  the  voltage 
of  the  mean  SOC  was  used  together  with  the  cycle  test  temperature 
of  35  °C.  All  discussion  of  cycle  aging  in  this  chapter  refers  to  these 
adjusted  values.  A  selection  of  typical  ’pure’  cycle  aging  curves  can 
be  found  in  Fig.  10. 

Within  literature  sometimes  a  linear  aging  over  charge 
throughput  Qis  found  [16],  but  also  work  about  square  root  aging 
was  presented  [7,17],  In  our  tests  we  found  two  different  de¬ 
pendencies  for  capacity  and  resistance.  For  capacity  loss  a  clear 
square  root  function  can  be  found  while  the  resistance  increase  is 
linear.  This  indicates  that  during  cycling  two  different  aging  pro¬ 
cesses  occur  in  the  battery,  one  that  has  a  strong  impact  on  capacity 
and  another  that  influences  resistance. 
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=  i-/vVa 

(18) 

=  1  +  /WQ 

(19) 

Analysis  was  done  with  a  two  step  fitting  like  before.  In  the  first 
step,  all  single  tests  were  fitted  using  Equations  (18)  and  (19).  After 
that  the  dependency  of  cycle  aging  factor  /?  on  cycle  depth  and 
average  voltage  were  analyzed. 

4.1.  Dependency  on  cycle  depth 

Cycle  depth  is  an  important  point  to  battery  aging.  For  a  lot  of 
battery  types  and  chemistries,  deep  cycles  cause  disproportional  or 
exponential  stronger  aging  compared  to  small  cycle  depths  [18],  In 
these  cases  oversizing  of  battery  capacity  is  an  effective  but  also 
expensive  option.  With  a  good  knowledge  of  cycle  depth  de¬ 
pendency,  the  oversizing  can  be  reduced  to  a  minimum. 

Fig.  11  shows  the  aging  factors  0  from  the  first  step  fits  over  cycle 
depth  ADOD.  Tests  at  equal  average  voltages  are  plotted  in  the  same 
color.  The  plot  indicates  a  non-exponential  dependency,  instead  a 
more  linear  trend  can  be  found.  For  capacity  loss  this  linear  de¬ 
pendency  seems  to  be  true  for  the  full  range  while  the  resistance 
increase  shows  a  constant  minimal  aging  for  low  cycle  depths. 

4.2.  Dependency  on  average  voltage 

Fig.  10  already  indicates  that  cycling  around  a  middle  voltage 
leads  to  the  lowest  aging.  Both  for  capacity  loss  as  well  as  for 
resistance  increase,  the  cycling  between  45  and  55%  SOC  has  the 
smallest  effects.  Cycling  around  lower  or  higher  voltages  increases 
the  cycle  aging  rate. 

A  plot  of  all  cycle  aging  factors  (I  over  average  voltage  is  shown 
in  Fig.  12.  For  all  cycle  depth  a  U-shaped  function  can  be  found  with 
a  minimum  around  3.7  V  average  voltage. 


For  lifetime  optimized  charging  strategies  the  dependency  on 
average  voltage  has  to  be  taken  into  account.  From  calendar  aging 
aspects  the  lower  SOC  has  advantages,  but  with  cycling  the  opti¬ 
mum  might  shift  to  a  SOC  level  around  50%. 

4.3.  Total  fit  Junction 

With  the  aforementioned  linear  cycle  depth  (ADOD  in  the  range 
of  0-1 )  and  quadratic  average  voltage  ( 0V  in  Volt)  dependency  a  fit 
was  made  to  all  aging  factors  /?  from  the  cycle  aging  tests.  The 
resulting  functions  for  cycling  capacity  loss  and  resistance  increase 


(3cap  =  7.348  •  lO-3-(0V-3.667)2  +  7.600  10  4 

+  4.081  •  10-3  ■  ADOD  (20) 

and 

/3res  =  2.1531O-4(0V-3.725)2  -  1.521  10  5 

+  2.798  -10  4  ■  ADOD  (21) 

with  all  values  of  /3res  below  1.510  5  being  set  to  this  value  to 
obtain  a  minimal  aging  as  seen  in  Fig.  11. 

The  total  aging  function  is  then  the  superposition  of  calendar 
and  cycle  aging  as  shown  in  Equations  (22)  and  (23).  In  this  func¬ 
tion  the  time  t  must  be  used  in  days  and  the  charge  throughput  Q  in 
ampere-hours.  With  this  function  a  calculation  of  battery  aging  and 
lifetime  is  possible  for  all  conditions. 

Cel  -  Ocap-t0'75  -  Patp's/o.  (22) 

R  =  1  +«res'ta75  +  /WQ 


(23) 
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average  voltage  /V 

Fig.  12.  a)  Aging  factor  /3cap  and  b)  aging  factor  (lres  over  average  voltage  for  different 
cycle  depths.  The  OCV  at  mean  SOC  of  the  cycling  was  chosen  as  average  voltage.  All 
cycle  tests  were  done  with  1C  at  35  °C. 


5.  Model 

To  calculate  lifetime  and  battery  aging,  the  voltage  and  tem¬ 
perature  of  the  battery  are  needed  in  each  time  step.  Most  appli¬ 
cations  only  provide  a  current  or  power  profile.  Therefore  a  holistic 
model  was  built  which  consists  of  three  parts: 


Fig.  13.  Graphical  representation  of  the  holistic  model.  Calculation  is  done  in  a  loop 
with  load  profile,  impedance  parameters  and  aging  factors  as  input.  Output  is  the 
remaining  capacity  and  resistance  together  with  a  lifetime  estimation. 


From  this  data  cycle  depth  and  average  SOC  are  calculated.  The 
average  voltage  for  cycle  aging  is  only  calculated  during  current 
flow,  rest  periods  do  not  influence  it.  Calendar  and  cycle  life  aging  is 
calculated  separately  for  every  single  time  step  using  the  aging 
functions  gained  from  the  accelerated  aging  tests.  To  get  the  total 
aging,  calendar  and  cycle  aging  are  added. 

The  impedance  parameters  are  then  updated.  All  capacitances 
from  the  equivalent  electric  network  are  scaled  like  the  normalized 
capacity  and  all  resistances  are  scaled  like  the  normalized  resis¬ 
tance.  After  that  a  new  loop  is  started  with  the  aged  parameters.  A 
schematic  of  the  aging  calculation  is  given  in  Fig.  15. 

6.  Verification 


•  Impedance  based  electric  model 

•  Thermal  model 

•  Aging  model. 

These  three  model  parts  are  called  in  a  loop.  A  schematic  of  this 
integrated  model  is  shown  in  Fig.  13.  As  the  model  is  impedance 
based  it  does  a  fast  computation  and  can  easily  be  parameterized 
through  electrochemical  impedance  spectroscopy  (E1S)  measure¬ 
ments  [19,20],  The  aging  model  is  based  on  the  equations  devel¬ 
oped  from  the  accelerated  aging  tests.  Every  aging  parameter  was 
analyzed  separately  and  put  into  a  mathematical  equation,  there¬ 
fore  the  effects  can  be  extrapolated  to  ranges  that  have  not  been 
tested  (e.g.  lower  temperatures). 

The  model  takes  a  load  profile  (can  be  power  or  current)  as 
input.  Furthermore  a  profile  for  ambient  air  temperature  can  be 
included,  else  a  constant  temperature  is  chosen. 

The  impedance  based  electrical  model  uses  an  equivalent  elec¬ 
tric  network  that  consists  of  a  series  resistance,  two  ZARC  elements 
and  an  OCV  source.  The  network  can  be  found  in  Fig.  14.  The  single 
elements  are  parameterized  using  lookup  tables  for  different 
temperatures  and  SOC  obtained  from  fitted  EIS  measurements.  The 
model  simulates  the  voltage  and  SOC  for  any  given  load  profile. 

A  thermal  model  is  coupled  to  the  electric  model.  Heat  pro¬ 
duction  is  based  on  ohmic  losses  calculated  by  the  electric  model. 
Heat  transfer  to  the  ambient  air  is  modeled  as  heat  conductivity 
with  variable  factor  to  match  different  mounting  scenarios. 

Voltage,  temperature,  SOC  and  current  are  taken  from  the 
combined  electric-thermal  model  as  inputs  for  the  aging  model. 


Along  with  the  accelerated  aging  tests  a  verification  profile  has 
been  applied  to  the  batteries.  For  this  verification  profile  a  driving 
cycle  was  used  together  with  charge  and  rest  periods.  The  driving 
profile  was  taken  from  real  car  measurements  done  in  Aachen, 
including  city  and  highway  driving.  The  power  was  scaled  down  to 
meet  the  specification  of  the  battery.  An  overview  on  the  driving 
part  of  the  profiles  can  be  found  in  Fig.  16,  the  power  together  with 
the  SOC  and  voltage  for  a  new  cell  are  shown. 

Three  different  test  profiles  were  built  with  this  load  profile. 
Plots  are  displayed  in  Fig.  17. 

The  first  and  second  profile  consist  of  two  driving  parts  a  day 
separated  with  long  rest  periods.  The  first  profile  has  the  charge 
periods  just  after  the  driving  part  while  the  second  one  charges 
before  driving.  The  third  profile  is  an  alternation  of  driving  and 
charging  without  long  rest  periods.  The  charging  in  all  three  pro¬ 
files  is  performed  with  1C 

The  profiles  tests  were  performed  with  monthly  changing 
ambient  air  temperature.  The  temperatures  are  average  values  for 
Germany  from  2001  to  2010  [21  ].  A  diagram  of  the  temperatures  is 
shown  in  Fig.  18.  The  initial  month  for  the  verification  was  September. 


Rser  ZARC  1  ZARC  2  OCV 


Fig.  14.  Equivalent  electric  network  used  in  the  impedance  based  model. 
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Fig.  15.  Schematic  of  the  calculation  process  within  the  aging  model.  Input  values  are 
given  for  every  single  time  step,  calendar  and  cycle  life  aging  are  integrated  on  these 
single  aging  steps. 


Temperatures  inside  a  battery  pack  build  in  a  car  might  be 
higher  than  the  ambient  air  temperature.  Therefore  higher  tem¬ 
peratures  have  also  been  tested.  Tests  with  all  three  profiles  have 
also  been  performed  at  the  average  temperatures  plus  10  °C  and  at 
plus  20  °C. 

Checkups  on  these  cells  have  been  done  every  month  with  the 
same  method  like  the  other  cycle  life  tests  described  in  Chapter  2.3. 

The  measured  data  for  both  capacity  loss  and  resistance  increase 
for  the  verification  profiles  can  be  found  in  Fig.  19  together  with  the 
simulation  results. 

The  figure  shows  that  the  aging  predicted  by  the  model  matches 
the  measurements  well  for  capacity  loss  although  the  model  was 
only  parameterized  at  temperatures  much  higher  than  the  tem¬ 
peratures  applied  in  the  verification  profile.  The  extrapolation  using 
the  Arrhenius  Equation  (6)  seems  to  be  valid  down  to  temperatures 
at  0  °C  Also  the  assumption  that  there  is  no  or  only  a  small  tem¬ 
perature  dependency  for  the  cycle  aging  proves  true  in  this 
verification. 

In  comparison  the  results  for  profile  1  and  profile  2  are  similar 
which  is  surprising  as  they  have  completely  different  charging 
strategies.  These  strategies  lead  to  different  average  voltages  during 
the  rest  period.  For  a  new  cell  it  is  4.174  V  for  profile  1  versus  3.970  V 
for  profile  2,  both  at  +0  °C.  Using  function  ( 14)  this  is  a  factor  of  1.25 
in  calendar  aging  between  the  two  strategies.  But  as  calendar  aging 
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Fig.  17.  Three  power  profiles  over  time  for  model  evaluation.  All  are  based  on  the 
driving  cycle  shown  in  Fig.  16.  a)  Charge  after  driving,  b)  charge  before  driving  and  c) 
no  rest  periods. 


is  much  smaller  than  cycle  aging  for  most  temperatures  the  influ¬ 
ence  on  total  aging  is  low  and  both  strategies  look  similar. 

The  simulated  resistances  for  most  verification  profiles  are 
lower  than  the  measured  ones.  The  biggest  relative  deviation  is 
found  in  cases  where  the  simulated  resistance  increase  is  very 
small,  for  higher  aging  rates  measurement  and  simulation  match 
better. 

As  this  cell  shows  only  a  small  resistance  increase  in  comparison 
with  the  capacity  loss,  the  last  one  will  limit  the  lifetime  for  nearly 
all  application.  Therefore  the  correct  capacity  prediction  is  more 
valuable  than  the  resistance  simulation. 

The  reason  for  the  deviation  in  resistance  increase  might  be  high 
currents  used  in  the  driving  profile.  While  all  accelerated  aging 
tests  were  done  at  1C  the  verification  profile  had  currents  up  to  3C 
for  discharging.  This  might  increase  internal  stress  on  the  particles. 
Beside  the  filming  reaction  which  is  included  in  the  aging  model,  a 
second  effect  seems  to  occur  which  increases  the  resistance 
without  consuming  lithium.  Due  to  the  high  currents  there  is  a  fast 
volume  change  during  intercalation  and  de-intercalation.  A 
possible  explanation  for  this  second  effect  could  be  partial  loss  of 
interparticle  contact  during  the  volume  change. 

7.  Conclusion 

In  this  work  an  aging  model  for  an  18650  cell  with  NMC  cathode 
chemistry  was  established.  An  extensive  set  of  both  calendar  and 
cycle  aging  tests  has  been  performed  including  more  than  60  cells. 
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Fig.  16.  a)  Driving  part  of  the  power  profile,  b)  SOC  and  c)  voltage  plot  for  this  profile 
and  a  new  cell.  All  values  are  plotted  over  time.  The  profile  is  taken  from  real  car 
measurements,  including  city  and  highway  driving. 


Fig.  18.  Temperatures  used  for  model  verification.  For  each  month  an  ambient  tem¬ 
perature  is  given  to  the  simulation.  Values  are  average  temperatures  for  Germany  from 
2001-2010  [2lJ. 


J.  Schmalstieg  et  al.  /  Journal  of  Power  Sources  257  (2014)  325—334 


333 


O  capacity  measurement  O  resistance  measurement  - capacity  simulation  - resistance  simulation 


Varied  parameters  were  voltage/SOC  and  temperature  for  the  cal¬ 
endar  tests.  In  the  cycle  aging  tests  cycle  depths  and  average 
voltage  were  studied.  During  all  tests  checkups  with  actual  capacity 
and  inner  resistance  were  determined. 

Using  these  values  for  capacity  fade  and  resistance  increase 
mathematical  functions  have  been  fitted.  For  both  calendar  and 
cycle  aging  tests  a  two  step  fitting  was  done.  The  resulting  fitting 
functions  include  the  dependency  of  capacity  and  resistance 
on  time  respectively  charge  throughput  and  the  other  varied 
parameters. 

For  the  tested  cell  the  calendar  aging  could  not  be  modeled 
using  a  square  root  function  of  time.  The  data  indicated  the 
occurrence  of  a  second  main  aging  effect  besides  normal  SE1  growth 
which  has  a  linear  time  dependency.  As  both  effects  could  not  be 
separated  by  the  fit  a  t0,75  dependency  was  used.  For  temperature 
dependency  the  Arrhenius  equation  was  used.  On  this  cell  the  aging 
speed  increased  linearly  with  voltage  for  a  calendar  aging  test. 

The  cycle  aging  data  was  adjusted  by  the  calendar  aging  predic¬ 
tion.  Fits  were  done  with  two  different  functions  for  capacity  loss 
and  resistance  increase.  The  capacity  was  implemented  using  a 
square  root  function  of  charge  throughput  while  the  resistance  had  a 
linear  dependency.  Both  capacity  and  resistance  had  the  lowest 
aging  while  cycling  around  3.7  V  (50%  SOC),  lower  or  higher  cycle 
ranges  increased  aging.  The  effect  of  cycle  depth  on  aging  was  linear. 

The  fit  functions  were  implemented  in  an  aging  model  which 
was  coupled  to  an  impedance  based  electric-thermal  model.  Three 
different  verification  profiles  were  tested  on  the  cells,  each  at  three 
temperature  levels.  In  all  nine  cases  the  predicted  capacity  loss 
matched  the  measured  data.  For  low  resistance  increase  the  model 
underestimated  the  aging,  for  higher  resistance  increase  the 
congruence  is  quite  good. 

This  holistic  model  is  able  to  predict  the  aging  for  different  drive 
cycles  and  management  strategies.  It  can  be  used  with  current  and 


power  profiles  as  input,  also  ambient  air  temperature  profiles 
changing  over  the  seasons  can  be  included.  The  model  is  a  powerful 
tool  to  optimize  BMS  strategies  and  lifetime  for  all  applications 
using  this  cell. 
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